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ABSTRACT 


This investigation is an attempt to verify the results of a theoretical model, 
developed by Karg] and Lim, for the scattering of sound from a poro-elastic sphere 
embedded in a poro-elastic host. It is a follow-on to that conducted by LT. Theodore W. 
L. Huskey. Both monostatic and bistatic measurements were taken on two porous glass 
spheres composed of 100 and 500 um glass beads and on an aluminum sphere. The 
Poisson’s Ratio was calculated from the shear and Young’s moduli measured from a 
cylindrical sample composed of 300 um glass beads. This was used to calculate the bulk 
moduli for the porous spheres; the shear moduli had been previously measured by LT. 
Huskey. These and other material properties were used as input to the theoretical model 
developed by Karg! and Lim. 

The experimental data were compared to the theoretical values. Reasonably good 
agreement between the measured and predicted scattering was obtained for the aluminum 
and 100 um spheres. The measured scattering from the 500 {1m sphere was in poor 


agreement with the theoretical predictions. 
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I. INTRODUCTION 


A. BACKGROUND 


Because of the increased emphasis on littoral warfare and the threat of the use of 
buried mines in these areas, there is increased interest in the scattering of sound in fluid- 
saturated porous media. 

Maurice Biot (Biot, 1956a, b) developed a general theory for the propagation of 
elastic waves in a fluid-saturated porous media. Biot’s theory has been applied to the 
scattering of elastic waves from a saturated porous sphere in a saturated porous host 
(Kargl and Lim, 1993). 

In 1993, LT. Huskey performed experiments in an attempt to verify Kargl and 
Lim’s model. The results of his experiment were good when 10% frame damping was 


included in the calculations. Other calculations were less satisfactory. (Huskey, 1993) 


B. OBJECTIVES 


Kargl and Lim developed a numerical model to compute the scattering of sound 
from a saturated porous sphere. This research is an attempt to experimentally measure the 
scattering of sound from a saturated porous sphere and compare the results with the 


numerical values obtained from Kargl and Lim’s model (Kargl and Lim, 1993). 


C. EXPERIMENT OVERVIEW 


As stated above, the purpose of this research is to measure the scattering of sound 
from a fluid saturated porous sphere. The spheres which were used in this research were 
composed of borosilicate glass beads. The beads had a mean diameter of 100 um and 
500 jim. The beads were coated with a heat curing epoxy powder and then poured into 
cylindrical molds. These molds were then heated to cure the epoxy. The resulting 


cylinders were ground into spheres with a diameter of about 6.8 centimeters. Cylindrical 


rods composed of 100 [1m and 500 um glass beads were made at the same time as the 
spheres. These rods were used by LT. Huskey to measure the permeability, porosity, and 
shear modulus of the porous matrices. These values are summarized in Table 1 (Huskey, 


1993). 





Sample 100 um 500 um 


Permeability, k 6.53x10"? m? 5.7410! m2 


Porosity, P 0.309 0.321 
Shear Modulus, G 2.81x10° Pa 2.72x10° Pa 


Table 1. Material measurements taken by LT. Huskey. 















Scattering cidasuremenits taken for these spheres were compared to theoretical 
values determined by the model developed by Karel and Lim. This model requires eleven 
material parameters to determine the scattered amplitude. The scattering solid provides 
five of these material parameters: mass density, bulk and shear moduli of the solid, and 
bulk and shear moduli of the porous lattice. Three material properties are determined by 
the structure of the porous lattice: tortuosity, permeability, and structural factor. The final 
three properties are of the bulk fluid: mass density, and the bulk and shear moduli. (Kargl 
and Lim). The program provided by Kargl used the bulk fluid’s bulk modulus and 
viscosity to determine complex values for the bulk and shear moduli. 

In this experiment all the material properties were accurately known with the 
exception of the bulk moduli of the spheres. This could not be measured by LT. Huskey 
at the time of his research (Huskey, 1993) due to the small aspect ratio (length to 
diameter) of the cylindrical samples. For this experiment a new cylinder of 300 um 
diameter porous glass beads with greater aspect ratio was obtained so that the bulk and 
shear moduli could be determined. From these moduli the Poisson’s Ratio can be | 
determined which can then be used along with the shear moduli of the 100 and 500 um 


samples to determine their bulk moduli. Once these parameters are determined they can 











be used as input, along with the other known material properties, to Kargl and Lim’s 


theoretical model and the results compared with the experimentally measured values. 





Il. DETERMINATION OF ELASTIC MODULI 


A. RESONANT ACOUSTIC METHOD 


The elastic moduli of a cylindrical sample can be obtained using a resonant 
acoustic method (Garrett, 1990). This method employs a transducer bonded to each end 
of a cylindrical sample which should have a length to diameter ratio >>1. The transducers 
are made from coiled magnet wire and are attached with epoxy. These are used to set up 
flexural, torsional, and longitudinal standing waves in the samples. The bars are 
positioned so that the transducers attached to the ends are centered between the pole faces 
of strong magnets. Depending on the orientation of the transducers to the pole pieces, any 
of the three vibrational modes can be selectively excited. Figure 1, from (Garret, 1990), 
illustrates the position of the magnet faces for the torsional and longitudinal modes. The 
fundamental frequency of each mode can be determined based on the length and 
boundary conditions of the sample. These frequencies can then be used to determine the 
moduli of the sample. 

In the experiments conducted by LT. Huskey only the shear modulus could be 
measured due to the length to diameter ratio of the samples (about 3:1) which intensified 


electrical cross-talk between the attached transducers. 





Figure 1. Transducer placement for the torsional and longitudinal modes. 


B. ELASTIC MODULI MEASUREMENT » 


A sample made of 300 tm glass beads with a length of 19.7 centimeters and 
diameter of 2.6 centimeters (giving a length to diameter ratio of 7.6:1) was obtained and 
the resonant acoustic method was used to determine both the shear and bulk moduli of 
the sample. Approximately 2 meters of #32 magnet wire were attached to each end of the 
sample with epoxy resulting in 15 turns of wire in each transducer, Care was taken to 
minimize the amount of weight added by the epoxy. The sample weighed 153.52 grams 
and the transducers added 4.38 grams. The mass loading of the transducers was 
significant enough that an effective length of the bar had to be calculated using the 


following equations (Garret, 1990): 


Lg = x1 + = | (longitudinal mode) 


L 


P x1 + = | (torsional mode) 
M 


where L is the actual length of the sample, m is the transducer mass, and M is the mass of 
the sample. 

The sample was positioned so that each transducer was centered between the pole 
pieces of the magnets. One set of magnets were rotated 90° to minimize electrical cross- 
talk between the two transducers (one transducer was also offset 90° relative to the other 
when it was attached). A Hewlett-Packard 35665A Dynamic Signal Analyzer was used to 
drive one transducer through a Hewlett-Packard 467A Power Amplifier. The signal from 
the other transducer was then input to the Dynamic Signal Analyzer. Both signals were 
monitored on an oscilloscope to ensure that no distortion was occurring. The output from 
the Dynamic Signal Analyzer was then swept upward in frequency from 100 Hz to 20 
kHz. The resonances were noted and the frequency bracket around each one tightened to 


allow for more accurate measurement of the resonance. 























Modal Frequency Summary 
| Mode Number, Frequency (Hz), Normalized Frequency (Hz), 
n f, f/n 


Torsional [SCS SCS. 





7830 2610 
10260 2565 
Average 2567245 


De! 
3940 3940 
8080 4040 










11775 3925 
15985 3996 
| Average 3975353 


Table 2. Resonant frequency measurements. 


The first four modes of both the longitudinal and torsional modes were measured 
and averaged together. The results are shown in Table 2. The Young’s modulus (£) of the 


sample was determined using the following equation (from Garrett): 
E=4pL,,"(f/n) 
where p is the mass density, n is the mode number, and f.” is the frequency of the nth 
longitudinal mode. The shear modulus (G) was similarly calculated using (from Garrett) 
G=4pL,,?(f7 /n) 
where f,’ is the frequency of the nth torsional mode. 


The calculated Young’s modulus was 3.8310’ Pa and the shear modulus was 
1.6810’ Pa. These values are approximately 1/20 of the borosilicate glass used to make 


the glass beads. These were used to calculate the Poisson’s Ratio (Vv) by combining 


_ 3K-2G 
2(3K + G) 
and Kae 
3(1-2v) 
E 
to give Yo —-—- 
2G 










2.81x10° Pa 2.72x10° Pa 







Shear Modulus, G, 
(from Huskey) 


Bulk Modulus, K, 2.96x10° Pa 2.87x10° Pa 
assuming v= 0.14. 


Table 3. Material properties of the samples. 









where K is the bulk modulus. The resulting Poisson’s Ratio was 0.14. This is an 
unusually small Poisson’s Ratio and 1s possibly due to the epoxy used to bond the glass 
beads together. The epoxy may be stretching between the beads allowing the cylinder to 
elongate without much lateral constriction. 

Since the 300 um rod was made at a different time than the 100 um and 500 pm 
spheres, the shear moduli measured by Huskey and the calculated Poisson’s Ratio were 
used to calculate new bulk moduli for both spheres. These values are listed in Table 3 and 


were used as inputs to Kargl’s program. 











fit. BACKSCATTER MEASUREMENTS 


A. MEASUREMENT OBJECTIVE 


The objective of the backscatter measurement was to separate the scattered 
acoustic pressure from the multipath interference. Multipath interference was caused by 
surface reflections from the transmitter’s side lobes and from piping used for filtration of 


the tank’s water. 


B. EXPERIMENT SETUP 


The backscatter measurements were taken in a water-filled tank, in Spanagel Hall 
Room 025, measuring 7.3 meters in length, 1.6 meters in width and 2.0 meters in depth. 
The walls and bottom of the tank are covered with anechoic tile. Figure 2 shows the tank 
setup for the backscatter measurements. All components were aligned along the 
centerline of the tank as viewed from the top of the tank. 

A type F33 general-purpose directional transducer was used as the projector. This 


transducer is shown in Figure 3 from the USRD Transducer Catalog, April 1991. Its 


Type F-33 B&K 8103 Sphere 
Transducer Hydrophone 





Figure 2. Tank setup for scattering measurements. 
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Figure 3. Type F33 general-purpose transducer. 


design consists of a small inner array made of 12 PZT disks with a frequency range of 
15-150 kHz. The larger outer array consists of 64 PZT squares with a frequency range of 
1-50 kHz (USRD Transducer Catalog). The arrays can be used individually or wired in 
parallel. In this experiment the arrays were connected in parallel for added directionality. 
To ensure accurate measurements, the scatterer must be placed in the far field of 
the transmitter. The following equation was used to determine the limiting distance to the 


far field (r,,,,.): 


where d is the dimension of the transmitter and A is the wavelength of sound (Kinsler et 
al.). The lowest frequency used in this experiment was 30 kHz, giving a maximum 


distance to the far field of less than 0.25 meters. 
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The receiver used was a Bruel and Kjaer type 8103 hydrophone. Figure 4 shows 
an illustration of the receiver from its calibration chart. It has a frequency range of 0.1 to 
180 kHz. The receiver was suspended from the arm of a protractor-like device illustrated 
in Figure 5. The arm can be rotated in five degree increments and locked into place with a 
lock pin. The 180° position corresponds to sound being scattered from the sphere directly 
back at the transmitter. The receiver could be moved radially, in increments of 5 
centimeters, between 40 to 75 centimeters from the center of the protractor. A small 
clamp was attached to the cable of the hydrophone allowing the depth to be adjusted. The 


clamp rested in any one of several beveled holes in the protractor’s arm. A weight was 





Units are mm 


Figure 4. Bruel & Kjaer Type 8103 hydrophone. 


suspended from the bottom of the receiver to ensure that it hung straight down. The 
weight was suspended a few centimeters from the bottom of the tank. 

The target spheres were suspended in a fine net by a string. The string passed 
through a hole drilled in the bolt which was the pivot point of the protractor’s arm. This 
allowed the depth of the sphere to be adjusted to correspond to the depth of the center of 
the transmitter. The spheres were degassed to ensure that no air bubbles were trapped in 


the spheres during the experimental measurements. This was done by placing them in a 


11 





beaker of water and then placing the beaker under a bell jar. A vacuum pump was then 
used to evacuate the bell jar until the water began to boil. At this point the hose to the bell 
jar was clamped and the vacuum pump was turned off. Periodically the beaker and sphere 
were agitated to dislodge any bubbles adhering to the surface. The vacuum was held for 
approximately 24 hours to ensure no bubbles remained trapped in the spheres. Next the 
hose clamp to the bell jar was removed allowing the pressure to return to normal. The 
bell jar was carefully removed and the sphere transported to the water tank in the beaker 
full of water. The beaker was submerged in the water tank and the sphere removed thus 


keeping it submerged at all times. 








BASE MADE FROM 3/8" STEEL PLATE. 


HOLES BEVELED TO MINIMIZE 
RECEIVER MOVEMENT. 





ALL UNITS ARE IN CENTIMETERS 


Figure 5. Protractor device used to position receiver. 
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HP 3314A 
Signal Generator 


Nicolet Pro 30 
Digital Oscilloscope 


HP 467A Type F33 
Power Amplifier 


Transducer 


DSS5020 
Oscilloscope 








B&K 8103 
Hydrophone 


Ithaco 1201 
Pre-Amp 





Figure 6. Electronic equipment setup. 


Figure 6 is a block diagram of the electronics used for this experiment. A 1 Volt 
peak to peak sine wave was generated with the Hewlett-Packard 3314A Signal Generator 
which was then amplified to a 10 Volt peak to peak signal by a Hewlett-Packard 467A 
Power Amplifier. This signal was then applied to the Type F33 transducer. The signal 
could also be monitored by either the Nicolet Pro 30 Digital Oscilloscope or the 
DSS5020 Oscilloscope. The signal received by the Bruel & Kjaer 8103 hydrophone was 
input to an Ithaco 1201 Preamplifier. The output from the preamplifier was then analyzed 


by the Nicolet Pro 30 Oscilloscope. 


C. EXPERIMENTAL PROCEDURE 


The frequency range used for this experiment was 30 kHz to 150 kHz, 
corresponding to a range of ka of approximately 4 to 22, referred to the wave number in 


water. 


1. Monostatic Measurements 

For the monostatic measurements the sphere was hung below the protractor. The 
B&K hydrophone was positioned 75 centimeters from the sphere along the arm of the 
protractor which was positioned at 180°. The signal generator was setup to send out 30 
kHz bursts at a rate of five bursts per second. Each burst consisted of a 30-cycle sine 
wave. | 

The Nicolet Pro 30 was used to monitor the received signal. A trigger delay was 
used so that the trace on the oscilloscope would begin just before the arrival of the 
scattered signal. The delay was determined by the speed of sound and the distance 
between the transmitter, sphere and receiver. The sample rate was set so that the pulse 
would fill as much of the oscilloscope’s screen as possible without being cutoff. Careful 
attention was paid to ensure that the sample rate remained well above the Nyquist 
frequency of the received signal. Fifty bursts of the received signal were averaged 
together and then saved to floppy disk. Measurements were taken between 30 and 150 
kHz in 2 kHz steps. 

At this point the sphere was carefully removed from the tank to a bucket of water 
and the above procedure was repeated to obtain the multipath interference background in 
the tank. The next measurement to be taken was the incident signal on the sphere. Due to 
the design of the protractor, the B&K hydrophone could not be positioned in the same 
place as the sphere but had to be positioned 15 centimeters behind it. This displacement 
was taken into account in the calculations by using a 1/r signal fall-off in the far field of 
the transmitter (Kinsler et al.). The incident signal was measured at each frequency. 

Once the received signal and multipath interference background measurements 
had been taken, the saved waveforms were subtracted from each other, leaving only the 
scattered pulse. The subtraction process was performed using the Nicolet Pro 30 Digital 
Oscilloscope. Figure 7 shows the measured signal, multipath interference background, 


and resulting scattered signal from the subtraction process. 
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-0.0300-+ 
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2.400 2500 
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Time (ma) 


Scattered signal. 


2400 2500 2800 2700 2800 


Figure 7. Example of subtracting the measured signal from the multipath interference 
background. At top is the received signal, in the middle is the multipath interference, and bottom 
is the difference between the measured signal and multipath interference background. 
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Next a Fast Fourier Transform (FFT) was taken of the scattered and incident 
waveforms, using the Nicolet Pro 30, to determine the amplitude of the desired 
transmitted frequency component. The start and stop points of the FFT were chosen so 
that they were inside the start and end points of the transmitted pulse. This was done to 
avoid distortion caused by the HP467A Power Amplifier turning on and off, and to avoid 
any ring up and down of the Type F33 transmitter. Care was taken to ensure that a whole 
number of wavelengths were taken and that the start and stop points were as close as 
possible to a zero crossing to avoid leakage into adjacent frequency bins when the FFT 
was taken (Hewlett-Packard, Application Note 243, pp. 25-26). No windowing function 
was used with the above method. The results of these measurements are listed in 
Appendix A. The method of calculation of the tabulated results are discussed in Chapter 


IV Section A, Normalization of Experimental Data. 


2. Bistatic Measurements 

Bistatic measurements were conducted in a manner similar to the monostatic 
measurements. The difference was that instead of adjusting the frequency between 
measurements, the angular position of the receiver was adjusted. The receiver was started 
in the 90° position and then moved in 5° increments until the 270° position was reached. 
The sphere was then removed from the tank and the background measurements were 
taken. Next the receiver was moved to measure the incident signal level on the sphere. 
After these measurements were made the resulting waveforms were subtracted and FFTs 
taken to obtain the scattered and incident signal levels. These measurements were 
performed at the following frequencies: 30, 60, 90, 120, and 150 kHz. Results of the 
measurements are listed in Appendix B. 

For the aluminum sphere measurements were taken only at 30 kHz and 150 kHz 


since these were to be used as a gauge of the effectiveness of the experimental method. 
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IV. DATA ANALYSIS AND RESULTS 


A. NORMALIZATION OF EXPERIMENTAL DATA 


All of the experimental data was normalized to a distance of | meter from the 


sphere. This was done by using the following equation: 


. V xr 
Normalized Scattering = —s#ee4¢ 


incident 
where r is the distance between the sphere and the receiver and V. 


incident 
-y x R+0.15 


Viaties incident R 


is given by 





where V’,,..iden 18 the measured incident signal (approximately 15 cm behind where the 


sphere actually was) and R is the transmitter to sphere distance. 


B. THEORETICAL DATA 


The values for the theoretical data were obtained from two FORTRAN programs 
provided by Kargl. One program calculates the theoretical monostatic data results and the 
other the bistatic results for a saturated poro-elastic sphere in a saturated poro-elastic 
medium. 

The program requires almost 30 inputs. These inputs are the material properties of 
the external fluid and poro-elastic medium, both water in this experiment, and the internal 
fluid (water) and poro-elastic medium of the scatterer. The internal poro-elastic medium 
was either porous glass or the aluminum. Inputs to the programs can be found in 


Appendix C. 


Cc. COMPARISON OF THEORETICAL AND EXPERIMENTAL DATA 


1. Scattering from the Aluminum Sphere 

Figure 8 shows a comparison of the theoretical and experimental data for 
monostatic scattering from the aluminum sphere. Experimental data points are marked by 
>< and connected with a dashed line. The dashes are only to guide the eye and are not an 
attempt at a curve fit. After the experimental data were taken the sphere was weighed and 
measured to verify the properties for input into Karg!’s program. It was found that the 
aluminum sphere weighed approximately 30% more than it should. The sphere was not a 
solid aluminum sphere, it was an aluminum shell covering a core of some heavier 
unknown material. The main features in the experimental data agree very well with the 
theoretical data below 90 kHz except that the experimental values are shifted to the right 
by 5-10 kHz. Beyond 90 kHz the features match in relative position but the experimental 
data has a lower amplitude. Figure 9 shows the data adjusted so that it spans 20 to 150 
kHz. This shows a much better match between the theoretical and experimental values at 
frequencies below 90 kHz. Since the sphere was determined not to be solid aluminum 
these results are considered to be in reasonable agreement with the predicted values and 
show that the experimental procedure is valid. 

For the bistatic case, shown in Figure 10 and Figure 11, the structure and 
amplitudes do not agree well. Again, this is not unexpected since the sphere was not solid 
aluminum. However, because the lobing is clearly defined in the experimental! data, it 


was again considered that the experimental procedure was valid. 


18 












0.030 Experimenta) Amplitude 


—— Theoretical Amplitude 


0.025 


0.020 
wo 
aw] 
2 
= 
E 
< 
o 0.015 
oO 
N 
S 
E 
o 
z 
0.010 
x 
0.005 
0.000 
30.0 50.0 70.0 90.0 110.0 130.0 150.0 
Frequency (KHz) 


Figure 8. Normalized monostatic scattering amplitudes from the aluminum based sphere. 


0.0300 


- X- - Shifted Experimental Amplitude 


—— Theoretical Amplitude 
0.0250 





x | | 
0.0200 2 ’ ; rn | y 
— an x ] | i 
JAM \ ig 
| ! i 1 : ! 
< | | ee x | x ‘ 
3 0.0150 1] eee t | a ae 
|| has at iP | 
1) RRA Ae 
a0109 Fe oe (es 
‘| . x) t . 
7 i 
{ 


0.0050 


0.0000 


20 40 60 80 100 120 140 
Frequency (kHz) 


Figure 9. Normalized monostatic scattering amplitudes from the aluminum based sphere with 
data ‘stretched’ between 20 and 150 kHz. 
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Figure 11. Normalized bistatic scattering at 150 kHz from aluminum sphere. 
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Figure 12. Normalized monostatic scattering amplitudes from the 100 um sphere. 


2. Scattering from the 100 um Porous Glass Sphere 


Figure 12 shows the results for monostatic scattering from the 100 iim sphere. 
Most of the features in the experimental data can be found to correspond fairly well to 
features in the theoretical model up to about 80 kHz and less so above 80 kHz. Note that 
the amplitude of the experimental data is lower than the theoretical data. This agrees with 
the bistatic data, which consistently shows that the amplitudes of the main lobes are 
lower than the theoretical prediction even when the side lobes have similar amplitudes as 
shown in Figures 13 through 17. 
Figures 13 through 17 show the theoretical and experimental bistatic scattering 
data sets. There is good agreement for the 30, 60 and 90 kHz data sets except for the 


amplitudes of the main lobes. The 120 kHz data matches well and the 150 kHz data only 


Z| 





agrees in overall magnitude of the scattering. This 1s believed to be due to the increased 
sensitivity of the beam pattern to the material properties. In varying the values used as 
input for Kargl’s program, it was noted that the beam pattern varied slightly at 30 kHz 
and increased in variation as frequency was increased. The resulting beam pattern varied 


significantly at 150 kHz. This is the most likely cause of the disagreement. 





Figure 13. Normalized bistatic scattering at 30 kHz from 100 um sphere. 
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Figure 15. Normalized bistatic scattering at 90 kHz from 100 um sphere. | 





Figure 17. Normalized bistatic scattering at 150 kHz from 100 um sphere. 
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Figure 18. Normalized monostatic scattering from the 500 um sphere. 


3. Scattering from the 500 um Porous Glass Sphere 

Monostatic scattering for the 500 ym sphere is shown in Figure 18. It is possible 
to see some correspondence between features below 90 kHz but above this frequency no 
match can be found. 

Figures 19 through 23 show the bistatic results for the 500 um sphere. There is no 
correspondence between the experimental and theoretical values. It is believed that this 1s 
due to defects in the 500 um sphere. Prior to the experiment it was noted that this sphere 


has a 0.9 x 0.5 x 0.25 centimeter gouge in its surface. The surface is also less consistent 


in texture from one area to another. Also the surface is more susceptible to crumbling 
than the 100 um sphere. Finally, both spheres were weighed and measured and their 


densities calculated. Both were found to have a density of about 1550 kg/m’. This density 
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corresponds with the 100 um bar used by LT. Huskey for determining the material 
properties of the 100 tm sphere. However, the density of the 500 um bar used by LT. 
Huskey was 1513 kg/m’ (Huskey, 1993) which is about 2.5% lower than for the 500 um 
sphere. 

Based on the material condition of the 500 tum sphere and on the good results 
obtained with the 100 um sphere it is believed that the 500 um sphere is defective and 


that this caused the disagreement between the theoretical and experimental data. 





Figure 19. Normalized bistatic scattering at 30 kHz from 500 um sphere. 
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Figure 21. Normalized bistatic scattering at 90 kHz from 500 um sphere. 
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Figure 23. Normalized bistatic scattering at 150 kHz from 500 um sphere. 

















V. CONCLUSION AND RECOMMENDATIONS 


A. CONCLUSION 


The scattering of sound from an aluminum based sphere and from porous glass 
spheres composed of 100 zm and 500 lum diameter glass beads was measured. 
Monostatic and bistatic measurements were made for backscatter at 30, 60, 90, 120 and 
150 kHz. Experimental data were compared with theoretical data computed by a 
FORTRAN program written by Kargl. The program was based on the theory developed 
by Kargl and Lim for the scattering of sound from a fluid saturated poro-elastic sphere in 
a saturated poro-elastic medium. 

Material properties of the spheres were measured for input into Kargl’s program. 
The elastic moduli of a cylindrical bar composed of 300 zm diameter bonded glass beads 
were measured and used to determine the bar’s Poisson’s Ratio. This Poison’s ratio was 
used to determine the value of the bulk moduli of the 100 um and 500 um porous glass 
spheres. Other material properties of these spheres had been previously measured by LT. 
Huskey. 

Comparison of the experimental data to the theoretical from Karg!’s program 
yielded reasonable results for the aluminum based sphere. Very good results were 
obtained for the 100 im sphere. The measurements at the lower frequencies agreed more 
closely than those at the higher frequencies. This is due to the increased sensitivity of the 
spatial structure of the scattering at higher frequencies to slight variations in the material 
properties of the spheres. The main lobe was found to be lower in amplitude than 
predicted at all frequencies. No explanation of this could be found at this time. It is 
speculated that this may be the result of slight inhomogeneity in the composition of the 
sphere. At higher frequencies it may be due to error in the axial alignment of the source, 
sphere, and receiver. The 500 um sphere produced extremely poor results and it is © 


believed that this sphere is defective, specifically it is thought to be non-homogeneous. 
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B. RECOMMENDATIONS 


Scattering measurements at lower frequencies is recommended. This would 
require a larger tank than is present at the Naval Postgraduate School. This would allow 
measurements to be made where the fluid’s viscosity would affect the flow of the fluid 
through the pores in the spheres. 

For measurements at higher frequencies the material properties of the spheres 
need to be known to fairly high accuracy. For this to be accomplished new spheres should 
be made at the same as the time cylindrical samples are made. The cylindrical samples 
should have a length to diameter ratio greater than 15:1 to allow accurate measurement of 
the elastic moduli. The surfaces of the cylinders should also be machined so that they are 
consistent with the sphere’s surfaces. Also, for measurements at higher frequencies, a 
more accurate method for positioning the receiver is needed. Measurements should be 
taken at sufficiently small intervals to adequately define the structure of the side lobes. 

Finally, an investigation into the non-uniform variation of the scatterer’s material 
properties should be performed. This would determine their affect on the structure of the 
beam pattern from the sphere which may explain the poor results obtained with the 


500 Lim sphere. 
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APPENDIX A. MONOSTATIC DATA RESULTS 
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Scattered | Incident Scattered | Incident Scattered | Incident 
Frequency } Amplitude | Amplitude | Normalized # Amplitude | Amplitude | Normalized | Amplitude | Amplitude | Normalized 
(mV) (mV) Amplitude (mV) (mV) Amplitude (mV) | (mV) Amplitude 
a EY 
32 89.1 
34 
36 
38 
r 
42 697 | 272.7 | 0.0183 
44 
46 
8 | 560 | 2254 | 00179 
| 734 | 2609 | 0211 | 424 | 2589 | oo7 | 642 | 4202 | 0.0109 
52 
54 
| 10.78 | 3449 | 0.0228 | 
58 
| 956 | 3963 | 0.0173 
64 


13.83 9697 | 0.0206 
1741 | 5653 | 002 | 1099 | 5969 | ooise | 3008 | 1.0805 | 0.0199 
17.87 11847 | 0.0119 
23.18 13123 | 0.0131 
19.27 


70 
72 
74 





76 | 10.92 | a27.1_| 0.0004 
78 38.38 | 15786 | 0.0174 
82 286 | 11814 | 0.0017 
a4 
88 
92 
94 
100 
104 | 6.04 | 46056 | 0.0009 
10a | 80.72 | 6755.6 | 0.0101 
Eee 


a1 
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Corrected Corrected | Corrected | 
Scattered Incident Scattered Incident Scattered Incident | 
Frequency | Amplitude | Amplitude | Normalized }| Amplitude | Amplitude | Normalized f Amplitude | Amplitude | Normalized 
(kHz) (mV) (mV) Amplitude (mV} (mV} Amplitude (mV) (mV) Amplitude 
118.05 | 9,125.0 | 0.0093 
112.10 | 0130 | 0.0100 
123868 69008 | 00128 
122 
0.0163 
[192 | 7oee | 30324 | cores | — | 2000 | — | i941 | 58063 | 0.0147 | 
196 [| 4489 | a7ete | ore | 2030 | 26534 | 0.0055 | 12071 | 51420 | 0.0168 | 
4747.1 | 0.0185 
3,611.2 
146 | 2693 | 17965 | cots | i797 | 16022 | 0.0076 | 49.04 | 3,444 | 0.0105 
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APPENDIX B. BISTATIC DATA RESULTS 









































30 kHz 
Scattered Scattered Scattered 
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| 1463 | oois7 [ato | ooze [146 | 00142 | 
| 2051.68 | tote 156 | ost 140 | ong7 
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60 kHz 
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120 kHz 
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150 kHz 
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APPENDIX C. INPUTS TO KARGL’S PROGRAM 


Monostatic Inputs 
Aluminum 100 um 500 um 
External Fluid (Water) 
Density, p; 
Bulk Modulus, K, 
Viscosity, 7 
Internal Fluid 
Density, Pr 
Bulk Modulus, Ky 


998.665 998.665 
2.17293x10° 2.17293x10° 


998.665 
2.17293x10° 


2S 
Oo 
Oo 
ro 


2700 


8.078x10"° 
0.001 


998.665 
2.17293x10° 
0.001 


998.665 
2.17293x10° 


Viscosity, 7, 0.001 


0.001 
External Medium (Water) Po 
Density, p 998 665 
Solid Bulk Modulus, K, (2.17293x10°, 0.0) | (2.17293x10°, 0.0) 
Lattice Bulk Modulus, K, (2.25x105, 0.0) (2.17293x10°, 0.0) | (2.17293x105, 0.0) 


(1.0, 0.0) 
1.65 
0.999999 


(1.0, 0.0) 
1.65 
0.999999 


Shear Modulus, u 
Tortuosity, a 
Porosity, B 
Permeability, k, 
Internal Medium (Sphere) 


(1.0, 0.0) 
1.65 
0.999999 


— 
© 


— 
© 
ak 
© 


Density, p>, 2700 2231 2231 

Solid Bulk Modulus, K,, (8.078x10"°, 0.0) (3.5x10"°, 0.0) (3.5x10'°, 0.0) 
Lattice Bulk Modulus, Kg} (8.078x10'°, 0.0) (2.96x10°, 0.0) (2.87x10%, 0.0) 
Shear Modulus, Uy (2.677x10"°, 0.0) (2.81109, 0.0) (2.72x10°, 0.0) 


Tortuosity, 
Porosiy, 8 =d?SCC~t ae SSC*dCSC*~=‘t dC 
Miscellaneous) SSS 
Catto 
ae tt tot ot 
Cateq TS t0 10 800 
distanced ds 
hash SO 
ee 
Trend 
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Bistatic Inputs 
100 um | 500 um 


mi 
pate 
® 
= 
3 
= 
a 
=: 
ror 
= 
Py) 
mal 
@ 
— 


998.665 | 998.665 


Density, p; 998.665 











Viscosity, 1 0.0 | 0.0 0.0 
ai Ae a Sa a 
2700 998.665 998.665 


2.17293x10° 


Ww 
c 
x 
= 
9 
Q. 
= 
S 
gw 
A 
a 
o) 
~I 
fe.) 
x 
—_—. 
ro) 
5 
; 
ro) 
ron) 
= 


External Medium (Water) | 

 cemiy.n ONS |=. 508.665 

(2.17293x10°, 0.0) | (2.17293x10°, 0.0) 
Shear Modulus, 1 (1.0,00 | (10,00) | (1.0, 0.0) 
arg 
Permeability, k, WO OP 


Internal Medium (Sphere) 
Density, Po 


Solid Bulk Modulus, K (8.078x10"°, 0.0) (3.5x10'°, 0.0) (3.5x10"°, 0.0) 
Shear Modulus, {1 (2.677x10"°, 0.0) (2.6x10", 0.0) (2.6x10"°, 0.0) 


porosiy, fp oo 806 
Permeably, ke | 1x0 | ——esox10% 


vA 


[oa aE a 2 ST 
[ties | —~stoq 5 —| same 
Tareq «120000 30000 _ 
“radius | 0.0346 
ore A 


Experimental Kyo ro (2.96x10°, 0.0) (2.87x10°, 0.0) 
Experimental (1, (2.81108, 0.0) (2.72x10°, 0.0) 
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